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Abstract. We present calculations of nuclear structure effects to the Lamb shift in light 
muonic atoms. We adopt a modem ab-initio approach by combining state-of-the-art nu¬ 
clear potentials with the hyperspherical harmonics method. Our calculations are instru¬ 
mental to the determination of nuclear charge radii in the Lamb shift measurements, 
which will shed light on the proton radius puzzle. 


1 Introduction 

The proton radius extracted from recent high-precision measurements of the Lamb shift (2S-2P atomic 
transitions) in muonic hydrogen atom yuH [1,2] displays 7 standard deviations from results determined 
from eH spectroscopy or e-p elastic scattering. This discrepancy raises intriguing questions about 
lepton universality and challenges theories of standard model physics. To resolve the proton radius 
puzzle, the CREMA collaboration at Paul Scherrer Institute started new measurements of the Lamb 
shifts in other light muonic atoms {pD, //^He^ and ju"^He^) [3]. These measurements aim to extract the 
nuclear charge radii with extremely high accuracy, which will be benchmarked with current and future 
spectroscopy measurements in electronic atoms [4-6] and electron-nucleus scattering data [7, 8]. 

The accuracy of the determined nuclear charge radius depends not only on the experimental pre¬ 
cision of the Lamb shift, but also on the theoretical calculations. The theories relate the Lamb shift 
6ls in a muonic atom/ion to the charge radius of a nucleus by 


^LS - <5qED + ^OPE + ^TPE, (1) 

where bQEo is the contribution from the quantum-electrodynamics (QED) including the photon vac¬ 
uum polarization, muon self energy, and relativistic recoil corrections. The other two terms in Eq. (1) 
are corrections due to nuclear structure. The term linear to (r^^) is dominated by one-photon exchange 
between the muon and the nucleus (Fig. 1), where JTIqpe ~ mj{Za)'^ll2 with = mfjMxlim^ ■¥ Mx) 
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denoting the reduced mass in the muon and nucleus (Mx) system. dxPE represents the two 
photon-exchange (TPE) contribution (Fig. 1), which is a combination of elastic (dzem) and inelastic 
(<5poi) parts, dzem is proportional to the third electric Zemach moment (r^;,)( 2 ) by [9] 


bzem = -m4(Za)5(r^\)(2)/24 . (2) 

The inelastic part 6poi is the nuclear polarizability effect, where the nucleus in the muonic atom is 
virtually excited in the intermediate states by exchanging two photons with the muon. 



Figure 1. The muon-nucleus one-photon (a) and two-photon (b) exchange processes: The bulb in (a) indicates 
the insertion of nuclear electric form factor. The bulb in (b) denotes the excitation of nucleus in the intermediate 
states between two photons. 


The bottleneck in accurately extracting the nuclear charge radius in Lamb shift measurements is 
the theoretical uncertainty in determining bxpE, especially bpoi. To determine of ^He and ^He 
from measurements of their muonic ions to a 3 x 10“^ accuracy, bpoi, which strongly relies on our 
knowledge of nuclear structures and dynamics, has to be known within a 5% accuracy [3]. 

bpoi can be represented as a combination of energy-dependent sum rules of a series of nuclear 
response functions. One way to calculate bpoi is to use inputs from measurements of photo-absorption 
cross sections and quasielastic scattering data. However, the applications of this approach in pD [10] 
and yU^He^ [11-13] do not provide the accuracy desired by the Lamb shift measurements. The high- 
accuracy requirement can be achieved alternatively through calculations using state-of-the-art nuclear 
potentials. For example, the studies of nuclear structure effects in pD using Argonne uig (AV18) 
potentials and pionless effective field theory improved the accuracy to a ~ 1% level [14-16], which is 
one order of magnitude more accurate than previous methods. 

We calculated the nuclear structure effects in /rD using chiral effective field theory (jt'EFT) po¬ 
tentials [17]. Combining the modern inter-nucleon potentials with few-body methods, i.e., effective 
interaction hyperspherical harmonics expansion (EIHH) [18] and Lanczos sum rules (LSR) [19], we 
also calculated bxpE in [20, 21], and [22]. By comparing the results using different 

potential models, we show that the uncertainty in bxpE is dominated by nuclear physics. Our results 
are key to the determination of from future muonic atom measurements. 

2 Nuclear structure effects in yuD 

bxPE in yuD was calculated by Pachucki [14] with the implementation of AVI8 nuclear potentials [23]. 
The dominant contribution is the sum rule of the electric-dipole response function, with corrections 
from the higher-multipole response function sum rules, relativistic and Coulomb distortion effects, 
etc. This calculation achieved a ~ 1% accuracy by systematically analyzing neglected higher-order 
effects in atomic physics. Our recent work [17] provided an independent derivation of the formalism 
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for ^tpe, which is consistent with Ref. [14].^ The later work by Pachucki et al. [16] carefully analyzed 
the relativistic-recoil effects in dxPE, which indicate, overall, a few parts per mil correction. 

To fully analyze the uncertainty in nuclear physics, we used in Ref. [17] not only the AV18 poten¬ 
tial but also ;fEFT with a series of parameterizations [24, 25]. We studied the order-by-order chiral 
convergence and the regulation-cutoff dependence. Our study shows a 0.6% uncertainty from potential 
variations, which must be included in the total uncertainty of dxpE in yuD. With a further investigation, 
we include an additional correction results in Ref. [ 1 7] to be consistent with Refs. [15, 16]. 

This term is proportional to the proton Zemach moment, its contribution to 5 xpe in is related 

to the elastic TPE effects in muonic hydrogen dzemf/tH) by 

a5<r3)(2)/24 = femOuH) (3) 

The hadronic polarizability in /tX can be approximated from its effects in /tH by 

^ z\z + N) [m,0uX)/m,0u//)]^ . (4) 

However, this estimate is based on the hadronic polarizability effect involving a single nucleon, but 
does not include the interference with nucleon-nucleon interactions. A more accurate calculation of 
^pof^ in pD is done by Carlson et al. by using high-Q^ part of the deuteron scattering data [10]. The 
combination of and (5p“['''(juX) yields a total hadronic TPE contributions in /uD to be = 

-47(10) ke V. Through collaborative communications among independent groups, bxpE in /uD is agreed 
to be 1.709 + 0.020 meV, which is adopted in a recent review by the CREMA collaboration [26] and 
will be used for extracting the deuteron charge radius in their Lamb shift measurement. 


3 Nuclear structure effects in yu'^He^, /i^He+ and yt/^H 

To extend the calculations to other light muonic atoms, we applied the EIHH method to expand 
the nuclear Hamiltonians and wave functions in the hyperspherical harmonic basis. This method 
is suitable for Hamiltonians with local (AVI8) or non-local (yt'EET) potentials. The energy-dependent 
sum rules of the response functions are obtained by implementing the LSR method [19]. 

We calculated the nuclear structure effects to the Lamb shift in /t"^He^, including the inelastic part 
of bxPE due to nuclear excitations (5™^) [20] and the elastic part due to nuclear charge distributions 
(bzem) [21]. We have recently improved the treatment of nucleon size effects and we show in Table 1 
updated results of for fjD and } Eor the hadronic TPE effects, rigorous calculations using 

scattering data are not yet available. Therefore, we estimate such effects in //"^He^ using Eqs. (3) and 
(4). We find ^ -0.34 meV and » -0.54(12) meV with the error dominated 

by the uncertainty in neutron polarizability. bxPE in /uD and //"^He^ are summarized in Table 1. 

Table 1. Two-photon exchange contributions to the Lamb shift in light muonic atoms in units of meV. 



^Zem 

£nucl 

^pol 

atomic error 

£hadr 

^TPE 

bxPE 


-0.424(3) 

-1.245(5) 

+0.018 

-0.047(10) 

-1.716(20) 

/He+ 

-6.29(28) 

-2.36(10) 

±0.10 

-0.89(12) 

-9.54(34) 


* Small discrepancies between Refs. [14, 17] appear in the reduced-mass-dependent coefficients of higher-multipole contri¬ 
butions. These are resolved later by Pachucki et al in Ref. [16], which agrees with Ref. [17], 

-The np coixelated nucleon-size conections to TPE were calculated in [20] with an inaccurate assumption that proton and 
neutron behave the same except for Coulomb. With a more rigorous calculation, we provide updated results. 
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We have recently extended our studies of TPE effects to /i^He^ and These muonic atoms 
of mirror isotopes, which contain unequal neutrons and protons, may provide additional information 
for the puzzle by probing the differences between muon-proton and muon-neutron interactions. The 
calculations of ^tpe in jU^He^ and will be detailed in an upcoming paper [22]. 

4 Conclusion 

We illustrate ab-initio calculations of TPE contributions to the Lamb shift in light muonic atoms 
originated from nuclear structure. We provide the results for different muonic hydrogen and helium 
isotopes, which are important for determining the nuclear charge radii in Lamb shift measurements. 
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